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The interplay of pairing and other interactions is addressed in this work 
using a simple single-j model. We show that enhancements in pairing correla- 
^ I tions observed through studies of the spectra of deformed systems, moments of 

' inertia, changes in transitional multipole amplitudes, and direct calculations 

of the pairing component in the wave function, indicate that even without ex- 
Q , plicit matrix elements responsible for pairing, a paired state can still appear 

' from the kinematic coupling of pairing to deformation and from other geo- 

ly^ ■ metrical restrictions that are of extreme importance in mesoscopic systems. 

CN , Furthermore, we demonstrate that macroscopic transitions such as oblate to 

prolate shape changes can lead to strong dynamic enhancements of pairing 
^ ' correlations. In this work we emphasize that the pairing condensate has an 

■ important dynamic and kinematic effect on other residual interactions. 

o 
o 

o 

-5 : I. INTRODUCTION 



The fact that a large number of nuclei in their ground state or in their low lying excited 
states are paired is supported by an overwhelming amount of experimental evidence. This 

><" 

' includes observation of odd-even mass difference, appearance of a gap in the spectra, reduc- 
ed ! 

tion of the moment of inertia, and analog of the Josephson effect in pair transfer reactions. 
By pairing in this work we imply the attractive interaction between pairs of nucleons on 
time-conjugate orbitals. It is widely accepted that the pairing interaction is responsible for 
creating a superconducting paired state, however the realistic interaction is much more di- 
verse than bare pairing. The complex interplay of all interactions that still leads to a paired 
state is far from being understood. Coherence between the pairs or even larger groups of 
nucleons can be formed in different quantum states, furthermore coherence may appear in 
the particle-hole (p — h) channel with other components of interaction contributing to col- 
lective excitations (in nuclei - shape vibrations and giant resonances) and deformation of 
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the mean field. All these effects are expected to dynamically and/or kinematically effect the 
paired state. There are also incoherent components of the interactions that introduce the 
stochastization of dynamics, but still can be influenced by the presence of collective features 
in dynamics. 

The appearance of a paired state is traditionally attributed to the strong short-range 
residual interactions between nucleons. However, in realistic nuclear systems all interaction 
matrix elements are correlated. There is no pure pairing interaction. It is well known in 
the theory of superconductivity, and in application to nuclear structure it has been shown 
by Belyaev [|l| long time ago, that an interaction with only pairing matrix elements would 
contradict the fundamental principle of gauge invariance. Recent studies of systems with 
two-body random interactions indicate that survival of collective phenomena such as 
pairing in realistic systems may hinge on these correlations between different types of matrix 
elements. These studies have shown that the paired ground state does not appear in sys- 
tems governed by two-body random interactions; furthermore, even weak but uncorrelated 
interactions of non-pairing type are very destructive with respect to the pairing state. 

Mesoscopic nature is yet another important property of nuclear systems. It has strong 
effect on kinematics and geometry of collective modes, interplay between different excitations 
and phase transitions. Finiteness was argued to be one of the main reasons for existence of 
the superconducting state in realistic nuclei 0]. 

In this work we show that observed pairing effects in nuclei do not result just from 
strong pairing matrix elements. Paired state appears from a very complex interplay of all 
residual interactions and their dynamic and kinematic behavior. Throughout this work we 
use a simple single j-level model with only one species of particle in order to discuss this 
interplay. This model provides strong kinematic constraints and the clearly pronounced role 
of antisymmetry requirements. Pairing problem can be solved exactly in a single j-level and 
treatment of all other interactions is substantially simpler. 

In Sec. ^ we introduce and discuss the kinematics of interactions in the single j-level 



model. The main results of this work are presented in Sec. |T|, where we investigate the 



dynamics of paired systems, and using a perturbative treatment of non-pairing interactions 
in the basis of paired states ||^ discuss the renormahzing effects of the pairing condensate 
on other residual interactions, consider the stabihty of the pairing condensate and evaluate 
the applicability of the pairing-based treatment. We emphasize that a single- j system is 
very kinematically constrained and for a number of independent choices of interactions the 
seniority, the number of unpaired particles, remains conserved. The interplay of pairing and 
quadrupole forces is studied in the "pairing plus quadrupole" (P+Q) model. We introduce an 
important concept of kinematic pairing as specific pairing effects that appear from kinematic 
restrictions present in a mesoscopic many-body system; they also contribute to dynamics 
of nucleon-phonon interaction With numerical studies we show that those effects, 

ignored in the standard P+Q model, result in a significant enhancement of pairing and 
directly influence the observable quantities, such as energy spectra, moments of inertia, 
and intensities of multipole transitions. For the systems with a nearly half-occupied shell, 
where the transition from oblate to prolate deformation takes place, pairing can be further 
enhanced due to the fact that in average the spherical shape tends to be restored in this 
region. 



II. KINEMATICS OF RESIDUAL INTERACTIONS 

The mean field is recognized as one of the most effective approaches in study of quantum 
many-body systems. Along with the shape and symmetry properties of the average many- 
body potential, the mean field also determines quantum numbers of elementary excitations, 
the quasiparticles. Low-lying states in the system as well as the response to external per- 
turbations can be understood in terms of the quasiparticles and their interactions, which in 
the lowest order are just pairwise collisions, see Fig. |1|. 

Spherical symmetry of the mean field is present in many nuclei throughout the periodic 
table. With the use of a spherical basis we guarantee the exact angular momentum conser- 
vation, avoiding approximate projections. Although the further discussion can be presented 
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in a general form, we limit our consideration to a single j'-level, that is $7 = 2j ' + 1 -fold 
degenerate. The general rotationally invariant two-body interaction Hamiltonian in a single 
j- shell, 



(1) 



defines the scattering of nucleon pairs coupled to angular momentum L . 
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Rotational symmetry and Pauli antisymmetry here result in the limitation that L is even, 
otherwise the scattering amplitudes Vl can be arbitrary. The L = term is responsible for 
pairing, the interaction of pairs on time-conjugate orbitals, |1) = \j m) and |1) = (— — 
m) . The strength of pairing is determined by Vq . 

A. Interactions in the particle-hole channel 

The interaction in Eq. (0) was given in the particle-particle (p — p) channel. The 
interaction can also be presented in the particle-hole (p — h) channel. The nucleon hole can 
be defined via a canonical transformation C , 
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The multipoles, particle-hole pair states coupled to a particular angular momentum K , are 
defined as 



Mkk = E 



^ J K J ^ 



V 



-mi K, 7712 



with the property {A4kkV = {—)'^ -Mr -n ■ The lowest multipole operators with K = and 
K = 1 are related to the constants of motion, number of fermions N, and components of 
angular momentum operator , 
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where f2 = 2j + 1 . 

The algebra of pair operators on one level is given by the following equations 
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It follows from the last expression that the odd multipolarity multipole operators form a 
closed subalgebra. Another subalgebra relevant to pairing is formed by operators Pqo , Pqo 
and A^oo 5 it will be discussed in detail in the following section. 

The original Hamiltonian (|1]) can be expressed in terms of interacting multipoles: 



V, 



H = T.yLY.nAPLA = eN-Y.^Y.^K.MK.. 

LA X ^ K 

The transformation formulas between the p — p and p — h channels on one level are 
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Vk = {2K + 1)Y.{2L + 1) l^^^ )Vl 
L \ 3 3 K 



(8) 



K \ J J L 



(9) 



5 



(10) 



K 



This transformation, often attributed to Pandya p], was first seriously discussed from 
the viewpoint of underlying physics and practically used by Belyaev see also pi]| , P 
Schematically the transformation from the p — p to p — h channel is shown in Fig. |I|. 
4 1 
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FIG. 1. Schematic diagram showing two-body scattering process. Two channels, particle-particle 
(s-channel) and particle-hole (t-channel), are indicated by the horizontal arrow with the matrix 
element Vl and by the vertical arrow with the amplitude Vk , respectively, where L and K are 
total angular momenta in corresponding channels. 

Fermionic antisymmetry requires that pairs of fermions on one level couple to even an- 
gular momentum, therefore interaction (|I|) is defined by j ' + 1/2 independent parameters Vl 
with L = 0, 2, . . . 2j — 1 . This fact is obscured in the Hamiltonian in the p — h channel, where 
particle and hole can couple to any angular momentum. The number of independent pa- 
rameters is still the same, however instead of a simple limitations Vl = for L = 1, 3, . . . 2j 
in the p — h channel the constraints for Vk are given by linear conditions 

Vk = {2K + 1) ^ (-)^+^' r ^ ^ \vk'. (11) 
K' [j J K j 

Analogous relations are known in the macroscopic Fermi-liquid theory. 

It is convenient to introduce a projection operator G that acting on Vl projects out only 
its physical component, 
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evL = —f^VL. (12) 

A similar operator also exists in the space of Vk , which is just a linearly transformed set of 
interaction parameters, see Eq. (^. However, it is no longer diagonal: 

QVK = lvK + l{2K+l)Y.i-f^^' \ ^ ^ ^' \vk>. (13) 

^' [3 j K j 

The condition that all Vl with odd L vanish is equivalent to QVl = Vl , similarly in the 
p — h channel Eq. (|13D is equivalent to QVk = Vk ■ 

Eq. (pUD can be viewed as an eigenvalue equation, where the kernel can be brought to 
a symmetric form by a simple rescaling of Vk by {2K + 1)^/^. This eigenvalue problem can 
be resolved by separating the 6-j symbol in the kernel with the recoupling technique. This, 
though, does not lead to anything new, because as a result one obtains that among 2j + 1 
eigenvalues there are j + 1/2 eigenvalues that are zero and the same number of eigenvalues 
equal to one, which is a consequence of this operation being a projection. Any physical 
eigenmode for the set of V corresponds to one particular Vl, and can be constructed using 
Eq. (W), since the projection operator is diagonal in the p — p channel. 



The special cases of Eqs. (fiA] ) and (|]) result in 



Y.Vk = E(2^ + ^)Vl = -nVo , E = -nVo , (14) 

K L K 

and effective single-particle energy in ([1^) can be written as 

These constraints are usually not addressed in nuclear models, because the Hamiltonian 
given in the p — h form is still good even if they are not satisfied, it merely contains the 
components that identically vanish in any fermionic many-body state, and still only j ' + 1/2 
independent combinations of parameters determine the interaction. The arbitrary amount 
of these components make the p — h form of the Hamiltonian expressing the same interaction 
not unique. The unique form can be reached if all non-physical components are removed 



with projection operators (|T^) or ([T3|). After projection the interaction remains physically 
identical to the original one, but the new parameters satisfy Eq. (0). 

The situation can be illustrated by an example of monopole interaction, where all nucleon 
pairs interact with identical strength, Vl = 1 for all even L , 

H= E EpIaPla. (16) 

L=0,2,...2j-1 A 

This interaction is very simple because its effect is only in counting the number of particle 
pairs in the system, therefore all states have the same energy 

^= \ ■ (17) 
Going to the j) — h channel Eq. (^, we rewrite (|16D in the form of interacting multipoles 

ii = ^N -Y7-^Y.^kMk. (18) 



where 



if 2 ^ 



Vk = -- [oko ^ — ) , and e = (19) 



All components of this interaction respect Pauli principle and Eq. ( [Tl| ) is fulfilled. However, 
it is not obvious that the action of this Hamiltonian is equivalent to counting the particle 
pairs. In order to gain a simpler form we add to this Hamiltonian a non-acting part 

H' = H+ E E^Ia^^a- (20) 

L=l,3,...2j A 

Transforming H' , where all = 1 , L = 0, 1, 2, . . . 2j , to the p — h channel, we get 

Vk = -^Sko, and e = ~. (21) 

Thus, only the monopole term is present in this interaction and 

H'.^. (22) 

Although Hamiltonians H and H' have a very different form in the p — h channel, they are 
identical in their action on a physical state. Despite the fact that introduction of inactive 
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components may allow for a simpler form of the Hamiltonian, the form where Eq. (|TT]) is 
satisfied, is preferred, not only because it allows to define interaction in the unique way, but 
also because it explicitly shows couplings between different physical p — h excitations by 
virtue of Eq. (pi]). 

The role that each interaction parameter Vl or Vk plays in determining the state of a 
many-body system is very complex, and generally for realistic systems these parameters are 
correlated by their common physical origin (such as core polarization or meson exchange for 
example) beyond the previously discussed kinematic restrictions. In realistic systems there 
are some Vl, Vk , and possibly their certain linear combinations, that have a significant 
tendency to form nuclear states with special coherent properties and symmetries. The Vq , 
pairing matrix element Vq, is known to be responsible for collective and macroscopic coherent 
effects similar to superconductivity and superfluidity in large many-body fermionic systems. 
Similarly, in the particle-hole channel V2 plays an important role for formation of collective 
vibrations and quadrupole deformation. 



Residual interactions can be formulated in the hole-hole (h-h) channel. With the trans- 
formation in Eq. (^) pair operators transform as 



CPLC-'=PL.^{-y-\ and CA^LC'"' = -ML, Ky^O. (23) 



B. Particle-hole symmetry 



The number operator transforms as 



CNC 



-1 



n-N. 



(24) 



With the help of the identity 




A 



(25) 



the Hamiltonian in Eq. (|I|) can be transformed to the hole-hole representation 
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(26) 



L A 



Since the number of particles (or holes) is a constant of motion, the p — p to h — h transfor- 
mation simply results in a constant shift of energy, while leaving the interaction invariant. 
The same result can be traced using the multipole-multipole {p — h) representation of in- 
teractions. Here all multipole terms with K ^ are invariant, and any changes are due to 
the monopole and single-particle terms. Particle-hole invariance results in important conse- 
quences [|T2|: an expectation value of any odd multipole moment of the A^-particle system 
is equal to the multipole moment of the corresponding state in the Q — N system; any even 
multipole moment is equal in magnitude but has an opposite sign in the corresponding states 
of the p — h conjugate system. In particular, the particle-hole symmetry requires that ex- 
pectation values of all even multipole moments identically vanish in the half-occupied shell. 
Therefore, a half-occupied shell can not be deformed. As we further show, this effect turns 
out to be helpful for preserving a paired state. This kinematic suppression of deformation 
is a result of a phase transition on the mean-field level from oblate to prolate deformation. 
Similar to the pairing phase transition, the mesoscopic nature of the system smoothens the 
sharp changes, thus extending the region of large fluctuations and suppressed deformations. 



The first steps towards understanding the nucleon pairing were taken even before 
Bardeen, Cooper and Schrieffer developed their powerful BCS method in 1957. The 
degenerate model involves a single degenerate single-particle level. The algebraic properties 
involving P = Pqo , -P^ = Poo 5 = TWoo operators on one level are particularly simple: 



III. PAIRING AND OTHER INTERACTIONS 



A. Pairing interactions and degenerate model 




(27) 



for even K , and 
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(28) 



for odd K , 

The important L = case, 



^ ^ 2v/2r+l 

^LO ^AO MlK 
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(29) 



(30) 



shows that the zero spin set of operators (P, P^ , and M) form the SU(2) algebra. By 
defining a quasispin £ , 



^ 2 4' 



'51 p, 

2 



(31) 



we can satisfy the commutation relations. The pure pairing interaction preserves quasispin; 
this can be converted into conservation of seniority, the number of unpaired particles s. This 
is the cornerstone of Kerman's |jl5| method and the EP algorithm |^ of exact solution of 
the pairing problem for the realistic level scheme. The eigenvalues of Cz and = C{C + 1) 
are related to the particle number N and seniority s according to 



N VL ^ n s 
r = r = 

^ 2 4' 4 2 



(32) 



By repeatedly commuting pair operators we get 



n-i / 2r2-2 2N 
n\l 
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(35) 
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The last expression results from [[Mkk, P^], P^] = . 

Collective paired states (the condensate) can be built on any state \N = s, s) with s 
unpaired particles by the simple action of the pair creation operators (-P^) \Si s), resulting 
in the state with n pairs in a condensate \N = 2n + s, s) . The normalization of such a 
state can be obtained using the momentum algebra, or iteratively with the help of the 
commutation relations, 

iN ^ .IP- (Pt)" 1^ (iV ^ .liV ^ .) . (36) 

The seniority formalism is useful because it takes all unpaired states |s, s) as a foundation 
upon which all other states are uniquely built by adding a paired condensate. The simplest 
lowest non-zero seniority states are the s = 1 state = 1, l(jm)) = a]m|0) and the 
s = 2 state (-Pl^oa)^|0) = |2, s = 2 (LA)), both of which are normalized to unity with our 
definitions. 



B. Pairing-based treatment on non-pairing residual interactions 

In this subsection we will assume that the system is paired, the ground state has seniority 
s = (assuming even A^), and the lowest excited state has s = 2 . Using the paired states we 
will evaluate the contribution of all residual interactions to the energy, the EP+monopole 
method using the states with s = 2 we will discuss the behavior of two unpaired nucleons 
in the presence of the — 2-particle condensate, and address the validity of the initial 
assumption that in the presence of all residual interactions the system is still paired. 

In the lowest order of perturbation theory for the s = state we have to examine the 
expectation values of all terms in the Hamiltonian of Eq. (^ for the paired state | A^, s = 0) . 
Following the commutation relations ( pSj) we obtain 

E(A^,0|>lL-Mx«|iV,0) = ^^|^f|^, ir = 2,4,..., (37) 

where only even multipoles contribute. The K = case is proportional to the square of the 
particle number 
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{N,Q\M^M\N,Q) = —. (38) 



Using (|3^) it can be shown that 



(AT, OlP^^Pi^liV, 0) = ^!^-^, L^O. (39) 

The L = corresponds to the solution of pairing in the degenerate model, 

N - s 

(AT, s\P^P\N, s) = (n-N-s + 2), (40) 

where the number of particles in the state is = 2n + s . This expression is valid for any 
seniority s . The expectation value of the Hamiltonian in the paired state is thus 

Vp Nin-N) NjN - 2) 

-y n-2 n-2 ■ ^^^^ 

The same result can be obtained in the multipole channel. The result ( PD is of particular 
interest since here the exact expectation value of the full Hamiltonian on the paired wave 
function is just the sum of the pairing and monopole contributions. The treatment of energy 
within the "exact pairing plus monopole" (EP+Mq) approximation is therefore the lowest 
order perturbation treatment [Q. It is also important to mention that s = 4 is the lowest 
seniority mixed with s = by the non-pairing part of the Hamiltonian. This is related to 
the conservation of angular momentum. The state of two unpaired particles cannot have 
spin zero since two particles on a single j-level have only one spin zero state which belongs 
to seniority zero. 

To consider the states of higher seniority it is convenient to utilize the quasispin group 
properties. All states with a given seniority have the same expectation value of and the 
number of particles in the paired condensate is reflected only in the quasispin projection Cz ■ 
With the help of the SU(2) quasispin group all operators can be classified by their seniority 
selection rules, and any expectation value (A^, s|X|A^', s') can be related to a quasispin- 
reduced matrix element using the Wigner-Eckart theorem. These procedures are 
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discussed in detail by Talmi |jT6|]. From the previously discussed commutation relations it 
follows that for odd L , AiiA is a quasispin scalar, while the even-L pair operators -P^A' 
Pla , and Ml a can be combined in components of quasispin vectors, for L = they define 
quasispin via Eq. ( pTD The Hamiltonian is a mixture of scalar, vector and second rank 
tensor in quasispin space, 

H = no + ni + n2. 

With the aid of the multipole expansion, the components of the Hamiltonian can be explicitly 
extracted. Due to the symmetry properties, the product of two identical vectors cannot 
have a vector component, because the cross product of two equal vectors is identically zero. 
Therefore A^J^^TW^k with non-zero even values of K contain no quasivector component. 
Thus, the quasivector part is fully contained in the K = terms 

Hi = -Vo (n-^]= -2VoC, . (42) 



2, 

The quasispin-quadrupole parts, that are only present in terms AI^^A^^k with even K , 
can be separated by decomposing the product, for example 

{C'f = i [2{C'f + + C- £+) + i (A{C' f -C^C- - £+) . (43) 

scalar quadrupole 

Therefore 



y« + ^y'^^.r^^ (44) 



3 V4 



and the remaining part is a quasiscalar 



(45) 



The quasivector part is proportional to and can act only within a multiplet, generating no 
change in seniority. Therefore in all transitions generated by the Hamiltonian and leading 
to a change in seniority the quadrupole part 7^2 is the only active component, changing 
seniority by either 2 or 4 units. Using the Wigner-Eckart theorem for transitions generated 
by the second rank tensor in seniority, we obtain 

14 





H\N,s-A) 


(s,s 
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s, s - 4) 



{N - s + A){N - s + 2){n - N - s + 2){n - N - s + A) 



(46) 



{N,s 


H\N,s-2) n-2N 


{N - s + 2){n- N - s + 2) 


{s,s 


H 


s,s-2) Q-2s \ 


2{n-2s + 2) 



(47) 

The situation with the diagonal in seniority contribution is somewhat more difficult since 
both components, quasiscalar and second rank tensor in quasispin space, are active in this 
case 



{N, siO\H\N, siO) = {s, smH\s, sie)) - UN - s)%. 
6iN-s)in-N-s) 

{s, s{0\Ti-2\s, )} ; 



(48) 



{n-2s){n-2s -2) 

^ here denotes all other quantum numbers not related to quasispin which are needed to 
identify the state. It can be convenient to extract a quadrupole component using two states 
in the quasispin multiplet: a state with no paired particles and a state with the same 
seniority but with one condensate pair. It can than be shown [|1^ that 

n-2s 



^2 {{s + 2, s{0\H\s + 2, s{C)) - {s, s{0\H\s, s{C)) + 2\/o%') • 

(49) 



The previously obtained formula for the s = case, Eq. (|4l|), results from the following 
conditions 

Q{Vo + 2Vo) 



{N = 0,s = 0\H\N = 0, s = 0) = , {N = 0,s = 0\n2\N = 0, s = 0) 
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(50) 



The s = 1 expression follows directly from Eq. (^) supplemented with 

(1, lijm)\H\l, l(jm)) = , (1, l{jm)\n2\h l(jm)) = , ^-u; ^ ^^^^ 
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{N, lijm)\H\N, l(jm)) = ^— i- ((fi - iV - 1)^ - (N - l)Vo 



(52) 
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The answer here contains only the pairing and monopole terms. An extra particle influences 
the pairing condensate only through the Pauli blocking. 

Two unpaired particles above the pairing condensate behave differently, and their inter- 
action is strongly renormalized. For the s = 2 case we get 

{N,2{JM)\H\N,2{JM)) = Vj - {N - 2)Vo + (53) 

(^-2)("-^-2)f(„_8),l£ + V'„)-2V^,- 1 



(0-6)(n-4) p 2 ' 2J+ IJ ' 

This equation shows that unpaired particles interact in the channel with angular momentum 
L with a reduced strength 

because of the presence of the A^p-particle condensate. The reduction is proportional to the 
expectation value of P^P and has a parabolic dependence on A^p, resulting in the maximum 
weakening of the pair interaction by about a factor of 1/2 . Addition of two unpaired particles 
implies an extra blocking of pairing and therefore requires more energy as compared to the 
case of the two particles being paired. The additional energy comes from a two-quasiparticle 
excitation and is proportional to Vq . The nontrivial contribution from other interactions also 
enters through the monopole and Vj terms. 

To investigate the stability of the paired state we consider the separation energy of a 
particle pair from the condensate S = E{N, s = 0) — E{N, s = 2) . In the approximation of 
large fl and Ap , Eqs. ( |41| ) and (0) give 



The self-consistency of this treatment based on pairing requires that pairing be stable and 



S < 0. Eq. (pSD as a function of the condensate size Ap/fi has three extremum points. 
The two points at the edges of occupancy Np/Q = and Np/Q = 1 are equivalent due to 
particle-hole symmetry and result in the obvious condition 

Vo<Vl, for any L ^ . (56) 
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A non-trivial condition appears in the third point of extremum, for the half-filled shell 
Np/Q =1/2, where 

Vn < — ^ — for any even L ^ . (57) 

2 2L + 1 ^ ^ ^ 



Quenching of residual matrix elements in the p — p channel, according to Eq. ( |5^ is an 
important phenomenon, which can prevent non-pairing interactions, especially ones that are 
incoherent with respect to the mean field deformation, from destructing the paired state. 
However, as can be seen from Eq. (^), the multipole-multipole correlations can damage the 



pairing state, and the above pairing-based treatment may become inappropriate. In Sec. 
we will continue the discussion of interplay of pairing and coherent multipole modes. 

C. Seniority conservation and kinematics of interactions 

The one level model is very restrictive kinematically, and constraints are somewhat fa- 
voring pairing. Although the interaction on a single level is defined with the explicit use 
of J ' + 1/2 independent parameters, such as Vl with even L, there are only few indepen- 
dent linear combinations that result in a seniority mixing interaction. Besides an obvious 
pairing component Vq , it is also possible to find a number of interactions that produce no 
quadrupole part = . It follows from Eq. (|^) that this will happen if for an arbitrary 
even L 

V. + ^^O. (58) 

This condition results in linear equation 

Y.'^ll'Vl' = 0, where JCll' = Sll' + {2L' + 1) \ ^ ^ \ . (59) 

Since 

/C^ = 3/C , (60) 
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the kernel K, has only two different eigenvalues, 3 and 0. All eigenvectors corresponding to 
the zero eigenvalue of )C are independent solutions of Eq. (|58|). Solutions of Eq. (|58D can 
be obtained as 



Vr 



J j L 
m —m 



L = 0,2, 



2j-l 



(61) 



The corresponding parameters in the p — h channel are 



V, 



2K + 1 



K 
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2j 



(62) 



In the above equation the second term in the brackets identically vanishes for all even values 
of . It is clear that Eq. (HD is satisfied by (|l|) and (H). Not all solut ions generated 
by J ' + 1/2 different values ofm = l/2, 3/2, ...j are linearly independent, this in general 
allows for existence of some independent linear combinations of interaction parameters Vl 
that result in seniority mixing Hamiltonians. The number of linearly independent solutions 
of Eq. (pHD is2k + l — 6ro, where k and r are determined as j + 1/2 = 3k + r , where r=0,l,2 



is the residue. ( |6T| ) generates all these solutions, they correspond to the zero eigenvalue of 
the kernel /C , and result in Hamiltonians that preserve seniority. Furthermore, since the 
quasiscalar Tig and quasivector Tii parts of the interaction result in a trivial A^-dependence 
of the spectra as follows from (^81) , 

(63) 



EiN, s(0) = E{N = s, s(0) - Vo{N - s) , 



the relative spacings between states of the same quantum numbers including seniority are 
independent of N for interactions that satisfy Eq. (^). 

As a remark we note that the 5- interaction can be generated by Eq. ( |6T] ) with m = 1/2 
and thus conserves seniority [Q. In addition to all these {2k + 1 ~ 6ro) choices, there is one 
trivial and linearly independent (for j > 1/2) freedom of selecting Vq that also results in a 
seniority conserving Hamiltonian. To emphasize these kinematic limitations we just mention 



that for j < 7/2 any interaction conserves seniority and as it has been recently noticed in ||T7 
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for j = 19/2 , for example, out of 9 possible independent choices of parameters only 2 lead to 
seniority non-conservation. In the limit of large j only about one third of parameters result 
in seniority mixing Hamiltonians. Finally, for the linearly independents sets of interaction 
parameters corresponding to the eigenvalue 3, the resulting Hamiltonian is almost purely 
quadrupole in quasispin 

n, = H-Vo[^-Ny no = -^. (64) 

The fact that it is possible to find non-zero sets of interaction parameters that result in a 
vanishing quadrupole in quasispin component of the Hamiltonian is non-trivial, it is further- 
more interesting that the number of independent possibilities is large. This is a result of a 
very restrictive kinematics of interactions on a single j-level. 

For a random choice of the two-body interaction, even with all of the above kinematic 
constraints the probability of getting a seniority preserving Hamiltonian is negligible as well 
as the conditions of condensate stability ( [57| ) are not necessarily satisfied. Thus, although 
random systems exhibit trends similar to those encountered in realistic nuclei with pairing 
I8| , p^ , the numerical studies show no enhancement of pairing in the low- lying states of 



random Hamiltonians. It has been demonstrated that in the ground state wave function 
of a random Hamiltonian on a single j-level the s = pairing component appears on a 
statistical level, i.e. with the same probability as any other component allowed by symme- 
tries. Therefore it has been argued that the presence of regular pairing, a prominent part 
of realistic physics, is not reproduced in randomly interacting systems 0. The intrinsic 
feature of interactions describing realistic systems is the presence of correlations between 
different interaction parameters. These correlations, along with kinematic features, such as 
discussed above and other dynamic couplings, is what makes the pairing effects survive and 
even dominate in the low-lying states of many realistic nuclei. 



19 



IV. PAIRING PLUS QUADRUPOLE MODEL 



As discussed above, the most general Hamiltonian can be separated into three parts, 
quasiscalar, quasivector and a second rank tensor in quasispin. The perturbation theory 
based on pairing treats exactly all quasiscalar and quasivector components. Since all odd 
multipoles are quasiscalars, only those non-pairing interactions that can be expressed in 
terms of the multipole operators of an even order, starting from K = 2 , are of interest as 
the most orthogonal to pairing. This leads to the Hamiltonian 

H = GP^P- ^Y.M^kMkk. (65) 

The lowest possible K = 2 multipole is responsible for quadrupole deformations and is 
usually the most energetically favorable. Thus, we will further concentrate on the pairing 
plus quadrupole (P+Q) Hamiltonian as defined below, 

H = GP^P-^Y.^LM2k. (66) 

The physically relevant parameters correspond to attractive pairing G < and attraction in 
the quadrupole channel, X2 > . This Hamiltonian is interesting for a number of reasons: it 
accounts for both short- and long-range parts of the residual interaction of nucleons through 
pairing and quadrupole parts, respectively, and consists of two very different components. 
Each of them separately is known to be responsible for collective phenomena, however acting 
simultaneously they lead to an interesting interplay. The study of pairing versus deformation 
within P-l-Q model is usually carried out with Hartree-Bogoliubov (HB) technique In 
fact the model is often defined as an arena for application of the HB method ||2^, ignoring 
exchange terms and previously discussed kinematic limitations. Studying the P+Q model 
beyond the HB approximation will be our further goal. 



A. Kinematic pairing 



The interaction parameters of the Hamiltonian ( pBD do not satisfy Eq. ([IT]) and, as 
previously discussed, this Hamiltonian contains a non-physical part. As a result the fact 
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that quadrupole part contributes to pairing as well as to all other components in the p — p 
or p — h channel, and that pairing makes a contribution to the quadrupole part is not 
seen explicitly. In order to observe these kinematic couplings we will reduce the form of 
the Hamiltonian ( |66D with projection operators to a unique form where conditions (|Tl|) are 
satisfied. We rewrite the Hamiltonian ( p5| ) in the following form 

(67) 



K K 



L A 

removing the unphysical part with the projection operation, 

VK = QKX2-{-fi2K+l)^. 

Thus 



(68) 



In particular, the monopole part is 




G 



(69) 







X2 + 2G 

2n 



(70) 



and there is a renormalization of the quadrupole strength which in the large Q limit behaves 

as 



X2 , X2 + '2G 



(71) 



In the particle-particle channel we have 



Vl = SloG + X2 



J J 2 



X2_ 

2n 



(72) 



J J L 

Eq. (|72|) shows that Vq = G — X2/^ , which indicates that even a pure attractive quadrupole 
interaction {G = 0) generates attractive pairing with the strength of the order Vq ~ 
Furthermore, as follows from the properties of 6j-symbols, even in the case of G = pairing 
is still the most attractive residual two-body force Vq < VljLq ■ Therefore, in this model the 
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pure attractive quadrupole-quadrupole Hamiltonian has a paired ground state (s = and 
J = 0) for nearly magic configurations N = 2 and N = Q — 2 , this effect was also observed 
in other studies of nucleon-phonon interactions 0]. A typical behavior of Vl and Vk as a 
function of L and K , respectively, is shown in Fig. 0. 

The above expression indicates that if Q is very large then only Vq ~ G and V2 ~ X2/2 
do not scale as l/fi, which leads to the usual P+Q model |^ with effectively decoupled 
quadrupole and pairing modes. This is not surprising because kinematic pairing as well as 
other kinematic couplings have mesoscopic geometry of nuclear systems as a source. 

In the P+Q model the condition for stability of pairing in a nearly full or nearly empty 
shell, Eq. (|56|), is fulfilled even without any explicit pairing component G = 0, since 
as discussed above Vq < Vl^q . However, an instability with the origin in the quadrupole 
channel appears in the middle of the shell, where Eq. ([57| ) leads to the following inequality 
in the limit of large Q , 



10 



(73) 
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FIG. 2. Interaction parameters in the particle-particle channel (upper plot) and particle-hole 
channel (lower plot) for a pure quadrupole-quadrupole interaction, X2 = 1 and G = 0, in the 
model space j = 15/2 are plotted as a function of multipolarity. Notice different scales on the two 
panels. 

For the general case of Hamiltonian (|65D, the condition (^) becomes 

G + — - < , K is even . (74) 

From the properties of the 6j-symbols it follows that the kinematic pairing resulting from 
any attractive multipole-multipole (even multipolarity) interaction is always attractive and 
is the strongest two-body component Vq < Vl+q in the p — p channel. The above result also 
indicates that lower multipoles [K = 2 is the lowest one) are more likely to destroy pairing 
because of the suppression factor 1/(2/^ + 1) . 

For completeness we present an exact equation for the separation energy of a pair from 
the condensate which is the minus excitation energy of the first J+ state with seniority s = 2 
for the model defined by the Hamiltonian (|65|) , 



S = G^Bj. (75) 

Here Bj is independent of G (J 7^ 0) and given by equation 

Y , (iV-2)(fi-iV-2) r,^ ^^^^ , , Ixj 



NY 

{Vt-2N + 2) 



where 



2(1] -2) 



Y = ^ E XK = G-Vo, (76) 



and Vj are determined as 



, j j K 

Vj = 5j,G + Y.XK 

^ \3 3 J 
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With pairing as the only interaction, the lowest excited state in the system with seniority 
s = 2 is at two-quasiparticle excitation energy, which is G in this case. Other interactions 
can lower this energy by a constant Bj which is mainly effected by the term xj/(2J + 1) . 
Since we are dealing with the zeroth order perturbation theory, this gap between the ground 
state and a lowest excited state of a given spin behaves linearly with the pairing strength. 
For example, in the P+Q model (see Fig. |]) the state J = 2 is the lowest excited state in 
the paired region. 




FIG. 3. The spectrum of lowest even-spin states relative to the lowest spin zero state in the 
system of six particles on j = 15/2 orbital for the P+Q interaction as a function of G . Two dashed 
straight lines that correspond to perturbation theory E2 — Eq = —G — B2 and E4 — Eq = —G — -B4 
are plotted for comparison (in this model B2 = 0.031 and B4 = —0.074). The quadrupole strength 
is set at X2 = 1 • 

The pairing based description for the system of six particles on the j = 15/2 level becomes 
unstable for G > —0.1 (negative G corresponds to attraction in the pairing channel) as 
predicted by Eq. (|73D with X2 = this agrees well with the comparison presented in Fig. 
^. However, in the region where unperturbed paired states can not be used to describe 
the system, such as the extreme case of no explicit pairing, G = , the effects of kinematic 
pairing are still quite strong. In order to emphasize this, the numerical values of overlaps 
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between the ground states of P+Q systems with no exphcit pairing (G = , X2 = 1) and 
paired systems {G = —1 , X2 = 0) are shown in Fig. ^. These results indicate significant 
pairing, that greatly exceeds the statistical prediction relevant to random interactions 0, 
shown in the figure by dashed lines. 
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FIG. 4. One-level systems with j = 15/2 and j = 21/2 , upper and lower plots, respectively, are 
considered with pure quadrupole-quadrupole interaction, G = . The wave function of the spin 
zero ground state is overlapped with the spin zero ground state of a paired system (G = —1 and 
X2 = 0), i.e. with the seniority s = state, that is unique for a single j-level system. The square 
of this overlap is plotted in solid line as a function of a particle number. This result is compared 
to the statistical expectation of pairing from random interactions, shown by dashed curve. The 
statistical expectation is defined here as an inverse number of spin zero states in the system [2]. 

The enhancement of pairing in the middle of the shell, observed in both plots of Fig. ^, 
is related to another kinematic feature. As it was discussed, exactly in the middle of the 
shell, due to the particle-hole symmetry the deformation must disappear which allows for 
more pairing correlations. In the following subsection we will further discuss this issue. 
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The same result is seen from the upper plot of Fig. ^ which shows that the pure 
quadrupole-quadrupole interaction, the region of G = , has a ground state dominated 
by the s = pairing component. Only in the region where kinematic pairing is explicitly 
balanced (this point Vq = is shown by the dashed vertical hne) by the repulsive explicit 
pairing given via positive G , the preference to s = pairing component disappears, and all 
J = states have almost the same overlap with s = wave function. The lower plot of Fig. 
^ addresses the low-lying properties of the spectrum in the same region of parameters; it 
will be discussed in the following subsection. 
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FIG. 5. Properties of the system of six particles on j = 15/2 orbital with the P+Q interaction 
are studied as a function of the parameter G ; the quadrupole strength is set at X2 = 1 • The upper 
plot shows the overlap of all six J = eigenstates in this system with the s = pairing state, 
which is defined as a ground state wave function for the pairing interaction G = —1 and X2 = . 
Lower plot shows the ratio {Ej=4 — Ej=q)/{Ej=2 — Ej=q) as a function of G, here -E'j=4,2,o are the 
energies of the lowest states with spin 4, 2, and 0, respectively. 



B. Pairing and deformation 

In this subsection we will investigate the interplay of pairing and deformation in the 
P+Q model. Rotational bands in the low end of the spectrum, and the Alaga intensity rules 
will serve us as tools in this study. The lower part of Fig. |^ shows that an indication for a 
rotational band near the "no pairing" region, G = . Here, judging by the excitation energy, 
the lowest states with J = 0, 2 and 4 are forming a collective rotational band. As it is clear 
from the figure this region is very small. Thus, for the most part we will concentrate here 
on a pure quadrupole-quadrupole interaction 

^ = -yE-^L-M2., (77) 

which, as we have shown above, is still significantly influenced by the kinematic pairing, 
with the bandhead J = state dominated by the s = pairing component. 

We note that although single-j model is very limited, an exact rotational spectrum can 
still be formed using 

H^X.ZMLM..^-^jf^J\ (78) 

which is seniority conserving and, according to Eq. (P) can be defined in the p — p channel 
with the aid of the set of parameters 

Vl = L{L + 1)-2jU + 1)^\ ^ ^ ^ ) . (79) 

[3 3 1 
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This interaction remarkably well satisfies the Alaga intensity rules, see below, for reason 
ably large j and . In the model with random interactions 0,0, the contributions ( |79| 



determine the effective statistical moment of inertia. Similar to this example other odd- 
multipolarity multipoles can be used to create seniority conserving interactions, this is an- 
other way of finding solutions to Eq. (^Sp. 

The main hint for the presence of deformation comes from the mean field approximation 
(with exchange terms ignored by definition of the model this is in fact a Hartree approx- 
imation). For axially symmetric deformation, the average values of the multipole moments 
are 

^ 2(3m^-j(i + l)) .... 
AI20) = V , rim, (80) 

m -1){Q^ -2) 

in terms of the occupation numbers Um = (ctm ) in the intrinsic frame with the 2;-axis 
oriented along the symmetry axis, and 

{M2-2) = {M22)=0. 

The single-particle energies in the body-fixed frame can be obtained via the usual self- 
consistency requirement 

2(3m^-j(j + l)) . .... 

Cm = -X2 I (-M20) • (81) 

This problem can be solved exactly |2ll]. The minimum energy corresponds to the Fermi 



occupation of the lowest pairwise degenerate orbitals \m\ = 1/2,3/2, .. .{N — l)/2 for 
prolate or \m\ = j, j — 1, . . . j — {N — 2)/2 for oblate shapes. The corresponding quadrupole 
moment is then given as 

for prolate deformation and 
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for oblate deformation. The deformation energy, defined as 

-^'deformation = |(-^2o)| ! (84) 

shows that oblate deformation is preferred for N < VL/2 and prolate correspondingly for 
> VL/2. In the middle of the shell there is a phase transition from prolate to oblate 
deformation. In this region energies associated with prolate and oblate deformation become 
equal. The average quadrupole moments at this point do not vanish; they are opposite in 
sign for prolate and oblate deformation. Because of this phase transition, mesoscopic nature 
of the system, and all other kinematic terms ignored in the model, the true ground state is 
a superposition of oblate and prolate forms in the region of a half-occupied shell. At the 
exact point of half-occupancy the particle-hole symmetry requires (A^2o) = for the true 
ground state of the system. The oblate to prolate transition turns out to be advantageous 
to pairing, as can be seen from Fig. |[ The content of pairing, the s = component in 
the wave function, is slightly enhanced in the middle of the shell. This enhancement, being 
accompanied by strong effects of kinematic pairing near both limits = and N = VL , 
is largely responsible for the presence of pairing correlations throughout the entire region 
within a single shell. 

The moment of inertia in the cranking approximation, that will be of use in our further 
discussion, is given by the following expression 

X=-E7^^b^l'- (85) 

With a sharp Fermi surface only four terms in the sum, |m| = |m'| ± 1 will work. Utilizing 
Eqs. ( pID and (0) we obtain 

_ ^(si^ - \){si^ - ^) _ ^](^]^-l)(^]^-4) 

folate , Xoblate 6 ^2 (^^ " iV)^ ' ^ ^ 

The Elliot's SU(3) model has a clear rotational structure and can serve as a link for un- 
derstanding a macroscopic deformation and its microscopic description. In the quadrupole- 
quadrupole Hamiltonian under consideration the SU(3) algebra is broken in such a kinematic 
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way that boosts competing pairing. The commutator of the quadrupole operators, according 
to Eq. (0), consists of vector and octupole components 



+14 E' 



2 2 1 



( 2 2 1 \ 



3 3 3 



2 2 3 



3 3 3 



2 2 3 



K K —q 



K K —q 



Mu 



(87) 



The octupole component is what breaks the SU(3) algebra; ignoring this term, and similar 
to Eq. (I), introducing a quadrupole operator according to 



M 



2q 



we obtain a SU(3) algebra in the standard form 

[Jq , Jq'] = —V2Cl^y^g, Jq+qi 
[Qq , Jq'] = -V^C'2g,^g' Q q+q' 



(88) 



(89) 



(90) 



[Qq ) Qq'] — 3"\/T0 Cag^^q/ Jq+q' ■ 



(91) 



The quadrupole-quadrupole Hamiltonian can be expressed via the bilinear Casimir operator 
of this group 



C = Q-Q + 3J- J. 



(92) 



The expectation value of this invariant operator depends only on quantum numbers A and 
fi that label representations, see for example p3| . 



(C) = 4(A2 + /i^ + a + 3A + 3/i) 



For a given representation (A fi) , where X > fi , the angular momentum can be 

K,K + 1,K + 2,...K + X K ^0, 



(93) 



A, A - 2 ... (1 or 0) 



K = 0. 



(94) 
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where an integer K takes values = /i, /i — 2, . . . (1 or 0) . Assuming a deformed band based 
on the ground state i^' = we choose = and A = N{VL — N)/2 where A coincides with the 
maximum possible value of angular momentum in the system. Thus in the approximation of 
SU(3) symmetry the expectation value of the quadrupole-quadrupole Hamiltonian becomes 



{N{J)\H\N{J)) 



(95) 



51](l]2_i) 51](fi2_i) 
The SU(3) model results in exact rotational bands with the A^-independent moment of 



inertia 



X 



12X2 



(96) 



In Fig. P the moments of inertia obtained from Eqs. (|86|) and (|96D and from fitting 
the J = 0,2,4 states in the spectra obtained from exact diagonalization are compared. In 
a view of the previous discussion it is not surprising that observed moments of inertia are 
significantly lower than ones predicted by both of the considered models. This effect should 
be attributed to pairing, which is mainly of a kinematic origin. 
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FIG. 6. Moments of inertia obtained using different models and from fitting the low- lying 
states in the exact spectrum are compared as a function of a number of particles. The pure 
quadrupole-quadrupole interaction is used (x2 = 1 and G = 0) for the system on a single j = 15/2 
level. The solid line corresponds to the moment of inertia coming from fitting the exact spectrum, 
short-dashed line corresponds to the moment of inertia from the mean field treatment, and long 
dashed line is for the SU(3) model, the latter predicts an A^-independent value. 

Static deformation of a nucleus manifests itself via relations (the so-called Alaga inten- 
sity rules) between the diagonal expectation values of multipole moments and off-diagonal, 
transitional amplitudes corresponding to the same K band. We consider here a quadrupole 
moment of the lowest excited J = 2 state and the E2 transition between this state and the 
J = ground state. In a single level model the quadrupole operator must be proportional 
to K as it is the only particle- hole operator with the correct rotational properties. We 
define quadrupole moment of the state as 

Q = {JM = J\M2o\J M = J) . (97) 

The B(E2) transition strength is defined as 

B(E2) = J2 \{JMf\M2.\JM,\\ (98) 

MfK 



The Alaga intensity rules are then expressed via relation p^j25 

Q2 4 



(99) 



B(E2) 49 

Within the pairing based treatment of interactions which implies no seniority mixing, 
this ratio can be calculated analytically. We assume that the lowest J = 2 state has pure 
seniority s = 2 . Due to the seniority conservation, this is the only state that can directly 
decay into the ground state, assumed here to have s = , via an E2 transition. The 
amplitude of this transition is given by 



(s = 2,AK|xLk = 0) 



^'^("-^> (100) 



\ {2K + i)n{n-2) 
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therefore in the hmit of strong pairing 

^^^^^ - n{n - 2) ■ 

The rate of this transition is maximized for a half-occupied system. The corresponding 
quadrupole moment in the J = 2 and s = 2 state can be determined from 

{N,s = 2,{JM)\Mko\N,s = 2,(JM)) (102) 



yM -M y 



J K J 



3 3 3 

Apphcation for quadrupole K = 2 results in the following expression 



2 = 1 (2£^£1£±1L. (103) 

7 Jn{D? -A) 



The quadrupole moment goes to zero for the half-occupied shell as required by the particle- 
hole symmetry. 

In general, the paired state is not deformed and thus the B(E2) transition probability 
in (|1U1D and the quadrupole moment ( |103| ) do not satisfy the intensity rule (^). With the 
quadrupole-quadrupole interaction presented, a deformation can appear and Alaga intensity 
rules can be fulfilled. In Fig. |^for the system j = 15/2 and = 6 , the quadrupole moment, 
B(E2) transition strength, and the ratio Q^/B(E2) are presented as a function of the pairing 
strength G . Dashed lines show the results for the strong pairing limit, obtained using Eqs. 
( |10T|) and (M. 
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FIG. 7. The quadrupole characteristics of the first excited 2+ state, the diagonal quadrupole 
moment Q, Eq. (97), middle panel, and the reduced transition probability to the ground state, 
B(E2), Eq. (98), upper panel, in the P+Q model for the single j = 15/2 level and = 6 , as a 
function of pairing strength G . Long dashed line on all plots indicates the results from pairing based 
description. The lower plot shows the prediction of the Alaga intensity rule (5^/B(E2) = 4/49, 
short dashed line. 

Fig. § compares the behavior of the ratio (5^/B(E2) as a function of the pairing strength 
G in different systems = 2,4,6,8, 10 with j = 21/2. There are two special cases. For 
N = 2 (the same is true for N = Q — 2) the ground state is paired, and the Alaga ratio is 
determined exactly via Eqs. ( |101D and ( |103| ). The second case is for = Q/2 , where due 



to particle-hole symmetry (5^/B(E2) = . For all other situations Alaga intensity rules are 
fulfilled to a sufficient degree of accuracy for the pure quadrupole-quadrupole Hamiltonian 
G = . As the pairing strength increases, the ratio (5^/B(E2) moves away from the Alaga 
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value towards the limit determined by pairing, which is shown by thin dashed lines. Curves 
corresponding to systems with a number of particles close toA'" = 2 (N — Q—2) or N = Q/2 
converge to the pairing limit significantly faster. This fact is yet another manifestation of 
preference for pairing correlations in these systems. 
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FIG. 8. Ratio (5^/B(E2) for systems with different numbers of particles on a single j = 21/2 
level is plotted in solid lines as a function of pairing strength G. Thin dashed lines, with cor- 
responding values of N marked on the right, indicate the intensity ratios for corresponding sys- 
tems in the limit of validity of pairing perturbation theory. All results should be compared with 
Q'^/B{E2) = 4/49 k, 0.082 , the Alaga intensity ratio of a rotating rigid body. 

V. CONCLUSIONS 

The pairing interaction is a very important part of the general residual interaction. 
However the fact that many nuclei are paired in their ground or low-lying exited states is 
a result of a non-trivial interplay of pairing matrix elements and other parts of the residual 
interaction, as well as kinematic constraints. It has been carUcr shown through a number 
of numerical studies that the realistic interaction is different from an arbitrary symmetry 
preserving random interaction. The difference lies in the correlations between the interaction 
parameters that exist even in the truncated shell model space, reflecting the overall properties 
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of the nuclear medium. Interactions that permit paired states and allow for deformations 
produce very correlated sets of residual two-body matrix elements. Apart from the dynamic 
correlations there are kinematic couplings between different collective modes, that mainly 
appear from the geometrical restrictions imposed on collective excitations. 

In many realistic systems pairing is relatively weak compared to the critical value of 
phase transition defined by BCS. However, as was shown by numerous authors, pairing 
correlations in mesoscopic systems extend far beyond the BCS phase transition. This makes 
the exact treatment of pairing a necessary component in the study of a sensitive interplay 
between pairing and other interactions. 

The simple single-level systems considered in this work has served as an excellent testing 
ground. The SU(2) quasispin algebra allows for an exact solution of the pairing problem, 
and perturbation theory based on the paired state p can be further developed with ease. 
This makes it possible to address important questions of the stability of the paired state as 
well as to investigate dynamic renormalizations of two-body interactions in the presence of 
the pairing condensate. 

The interaction parameters for a one-level system in the particle-particle and particle- 
hole channels can be related to each other, revealing kinematic constraints in a simple form. 
Additional restrictions on the behavior of the system come from the exact particle-hole 
symmetry. As a result of this symmetry, all multipole moments of even multipolarity are 
identically zero for a half-occupied shell. These constraints turned out to be very important 
for the preservation of pairing in the presence of deformation. Interplay of pairing and 
deformation was discussed using the P+Q model. The particular case of this model, the 
pure quadrupole-quadrupole interaction with no explicit pairing, has been addressed in 
detail. This interaction was expected to be the one most "orthogonal" to pairing; however, 
strong evidence of pairing correlations was found even in this case. We have shown that the 
pairing is the strongest attractive component in the quadrupole-quadrupole interaction when 
the latter is converted from the p — h to p — p channel. This means that for any attractive, 
i.e. favoring deformation, quadrupole-quadrupole interaction on a single j-level, the ground 
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state of a two particle or two-hole system is a paired state (J = and s = 0). The same 
is in fact true for a more general Hamiltonian containing any attractive multipole-multipole 
interaction of even multipolarity. This results in enhanced pairing correlations in near-magic 
configurations. An additional enhancement of pairing was observed in the systems close to 
the half-filled shell. This is related to the prolate-oblate phase transition taking place in this 
region; it weakens deformation thus allowing for more pairing. A number of calculations 
was performed to emphasize the discussed aspects, and the results were compared to the 
predictions of the mean field approximation and Elliot's SU(3) model. Direct overlaps with 
the s = paired wave-function, excitation spectrum, intensities of transitions and moments 
of inertia all indicate an appearance of pairing, consistent with the kinematic constraints. 

The simple model used in this work is only a very limited approximation to the situation 
in nuclear systems. However, the observed effects of pairing enhancement via kinematic or 
dynamic interplay with other interactions have to exist to some extent in realistic systems. 
The situation in the realistic shell model becomes much more diverse as it is complicated by 
an increasingly large number of interaction parameters, weakening of antisymmetry require- 
ments and other kinematic restrictions, and a diversity of collective modes. Pairing itself 
becomes different, isovector and isoscalar pairing modes can compete, the pairing state of 
the lowest total seniority s = is not unique and can involve coherent and incoherent pair 



vibrations p6[ which, under effects of other interactions, may result in a paired condensate 



that is different from the prediction of the regular BCS. 
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